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Abstract 

We study the surface shape of two-dimensional piles using experi- 
ments and a continuum theory for surface flows of granular materials 
(the BCRE equations). We first obtain an analytical solution to the 
BCRE equations with a simple transformation and show that the sur- 
face shapes thereby predicted are in good agreement with the experi- 
mental results. By means of such an analytical solution, we find that 
the formation of the curved tails at the bottom of such piles depends 
not only on the properties of the granular materials but also on the 
drift velocity of the grains within the rolling layer. 
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Granular materials have received increasing attention in recent years on 
account of the fascinating granular flow patterns and their important prac- 
tical applications. It is well-known that the angle of repose, a characteristic 
property of granular materials, is determined by measuring the steepest angle 
of the surface of granular piles. However, the surface of such piles, formed by 
pouring granular materials on a horizontal table, is not always the straight 
line implied by the constant angle of repose, in the sense that a curved tail 
is often present at the bottom of the piles. Alonso and Herrmann [TJ first 
investigated the shape of such a curved tail of a two-dimensional static sand- 
pile. By introducing the notion of kinks, they built a model to describe the 
local slope of the two-dimensional sandpile surface (cf. Eq. (1) in Ref. pQ), 
and then, using the observed translational invariance, derived the following 
expression for the shape of the tail, 



where 9 r = tan# c with 9 C being the angle of repose, h(x) is the height of the 
sandpile surface, x represents the corresponding horizontal coordinate, h m 
is the maximum height of the sandpile and l e {= l/r) is a fitting parameter 
with I being the horizontal length of the kink (which typically is of the size 
of a grain) and r the corresponding accumulating rate of the grains on the 
kinks. The above expression shows that the existence of tails introduces 
a logarithmic correction to the straight slope given by the angle of repose, 
and can be used to fit the experimental surface shape quite accurately by 




(1) 
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adjusting the magnitude of l e . However, the derivation of Eq. JTJ) entails 
three intermediate parameters, i.e. p, I and r (where p is the density of kinks 
at height h), but none of them can be determined experiment ally (cf. Ref. pQ 
for details), so that the final fitting parameter l e cannot clearly identify which 
factors have an effect on the formation of the curved tail at the bottom of 
the pile. 

In order to answer this query, we study the same problem described 
above on the basis of a continuum theory of surface flows of granular mate- 
rials proposed by Bouchaud, Cates, Ravi Prakash and Edwards (the BCRE 
model) J2). This model makes use of two coupled variables to describe the 
dynamics of the two-dimensional sandpile surface: the height of the static 
sandpile, h(x,t), and the thickness of the rolling grain layer, R(x,t), both 
of which are a function of the horizontal coordinate x and the time t. On 
assuming that the drift velocity, v, of the grains within the rolling layer is 
independent of x, and that diffusion effects are negligible, the original gov- 
erning equation of motion for the rolling grains in the BCRE model can be 
simplified into 

dR(x,t) = _ v dR(x,t) + 
ot ox 

where the interaction term V accounts for the conversion of immobile grains 
into rolling grains, and vice versa. Given that the deviation of the local 
surface slope 8(= —dh(x, t)/dx) of the pile from the critical value 6 r (= tan# c 
with 9 C being the angle of repose) is everywhere small, the simplest form of 
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r is, to a first approximation, 



T(h, R) = - 7 (0 r - 9)R{x, t) = -7 r + — ^ t), (3) 



where 7 is a positive constant relative to the property of granular materials, 
and can be interpreted as equal to the rate of collision between rolling and 
static grains. Moreover, the principle of mass conservation requires that 



In what follows, we shall first directly apply the simplified BCRE model, 
i.e. Eqs. (J2J) and (jUJ) (henceforth referred to as the BCRE equations), to 
investigate the surface shape of two-dimensional piles. 

As a granular pile reaches steady growth, the shape of the pile is found 
to be essentially independent of time (t), i.e. the front of the pile propagates 
along the x direction with a constant velocity (u). Therefore, we can combine 
the two independent variables, x and t, into one (say z) by means of the 
following transformation, 



dh(x, t) 
dt 



-r. 



z = x — ut. 



(4) 



Consequently, the BCRE Eqs. (J2J) and (j3J) become 





and 
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We can rearrange Eq. (jSJ) to yield 

dh(z) fv — u\dR(z) 



(7) 



dz \ u J dz 

which, when substituted into Eq. (jUJ), yields an ordinary differential equation 
for the thickness of the rolling grain layer R(z), 

- u) (l + ^R(z)^j = -j8 r R{z). (8) 

We integrate Eq. (|SJ) and, after applying the boundary condition which de- 
termines the flux of pouring grains (i.e. R = R m at z — 0), obtain the 
expression for R(z), 

z9 r = — In -§p- + l — (R m - R(z)) , (9) 
7 R{z) u 

where R m is the maximum thickness of the rolling layer. 

Similarly, on integrating Eq. (JZJ) subject to the the boundary condition 
{h — at R — 0), we find the relationship between R(z) and h(z), i.e. 

a 

R(z) = ——h(z). (10) 
v — u 

Substituting Eq. (fTUj) into Eq. Q, we therefore obtain the expression for 
h(z), 

h m -h(z) v-u h m 
where h m is the maximum height of the static pile which must satisfy 

h m = Rm- (12) 

U 
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Actually, the maximum thickness of the rolling layer (R m ) usually is much 
smaller than the maximum height of the static pile (h m ), which, in view of 
Eq. (JI2J), leads to 

v > 2u. 

Therefore, Eq. (fTTj) simplifies to 

h m -h{z) v h m 

r 9 r j h(z) 

which reduces to, 

h m -h(x) v h m 

x = a + T~ ln TT^> ( 14 ) 

u r Or'y n[x) 

provided that the growth of the pile is slow enough to satisfy the quasi-static 
assumption, i.e. ut ~ 0. 

It is interesting to note that Eq. (|14j) has the same form as Eq. (|1J) in 
that both equations have two terms on the right-hand side: the first one just 
corresponds to the linear part implied by the angle of repose 9 C , while the 
second one represents the non-linear part which predicts the existence of the 
logarithmic tail at the bottom of the pile. In fact, the only difference between 
Eqs. (0) and (fTlj) is the coefficient of the logarithmic term: one is l e and the 
other is vj^Q T . In contrast to the former (a fitting parameter in Eq. (JTJ), 
the latter consists of an operational parameter (v) and two phenomenological 
constants (8 r and 7) characterizing granular materials, all of which, in prin- 
ciple, can be obtained directly by experimental measurements [3J. Therefore, 
according to Eq. (JHJ), we clearly see that, for given h m , the length of the 
curved tail depends on the rolling velocity (v), the angle of repose (9 C ) and 



the rate of collision between rolling and static grains (7). Moreover, since v 
as well as 7 are relevant to the dynamics of the rolling grains on the top of 
the static pile, their presence in Eq. (J 14)) throws new light on the formation 
of the curved tail at the bottom of the pile. Thus, for a quasi-static pile, the 
flow of rolling grains is due to very small perturbations on the pile surface 
where the length, time, and velocity are respectively scaled with respect to a, 
I/7 and 07 (where a is the size of the grain) so that v must be proportional 
to 07. However, since 07 is essentially independent of the perturbations, so 
is the value of v. In this case, the formation of the tail is predicted to depend 
mainly on the size and the shape of granular materials. On the other hand, 
if the external disturbances are not so small, the value of v may be relevant 
to these disturbances rather than being simply scaled with 07. Under this 
condition, the formation of the tail depends not only on the properties of 
the granular materials but also on the operational condition, i.e. the drift 
velocity of the grains within the rolling layer. 

In order to check the reliability of Eq. (jl4|) . we also performed experiments 
with five different granular materials in a rectangular two-dimensional cell 
made of two vertical Plexiglas plates of size 22cm x 45cm with a fixed gap of 
2.5mm. Granular materials were poured through a funnel placed on the up- 
left corner of the cell. Several pictures of static piles are shown in Fig. and 
the qualitative comparison between them is listed in Table 1. For example, 
no obvious tail exists in the static pile of spherical glass beads (cf. Fig. 1(a) 
where a = 0.3mm and 9 C = 30°), but a long tail is present at the bottom 
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of the static pile of cubic sugar grains (cf. Fig. 1(c) where a = 0.8mm and 
6 C = 45°). These experimental observations of the static piles show that the 
formation of the tail is relevant to the properties of granular materials, which 
is consistent with the previous discussion following Eq. (JTJJ). 

As mentioned earlier, all the parameters in Eq. (jl4j) . i.e. v, 8 C , 7 and 
h m , can be obtained by experimental measurements j^j. However, compared 
with the last three parameters, it is more difficult to accurately measure v for 
the case of non-uniform particle velocity profiles within a varying thickness 
of rolling layer R(x,t). Under this condition, we can fit the experimental 
results by taking v to be the only adjustable parameter and using the values 
of 6 C , 7 and h m as determined from the precise experimental measurements. 
Here, 6 C and h m are measured in our experiments while 7 is obtained from 
early measurements jSj where the experimental setup as well as granular 
materials were similar to our case. In Fig. a quantitive comparison of 
the surface shape is made between experiment and theory: the dark part is 
the experimental sugar grain pile, and the solid line represents the surface 
shape given by Eq. JHJ). Clearly, the agreement is satisfactory in both linear 
and non-linear parts with v = 11.5cm and 7 = 23/s [3J. The value of v is 
consistent with previous estimates [Sj under the same flux condition as in the 
present experiments. 

Next, we investigate the dynamic behavior of pile growth using sand 
in Fig. 1(b). During the experiments with sand, we recorded the whole 
process of the sandpile growth, and thereby obtained a family of sandpile 



8 



surface profiles at different growth stages (cf. Fig. 0(a)). At the early growth 
stage (e.g. t = 13s), the corresponding sandpile was small and the tail 
was clearly seen; but when the sandpile grew larger (e.g. t = 37s), the 
corresponding tail became hardly discernible. Since the position of the funnel 
was fixed while pouring granular materials into the cell, the distance between 
the outlet of the funnel and the top of the sandpile, called H, decreased with 
the growth of the sandpile. Consequently, the initial velocity of the rolling 
particles also decreased with time, which may result in the decrease of the 
tail length. To examine this idea, we can still make use of Eq. (|14j) to fit 
the experimental profiles of the sandpile surface at different growth stages 
by choosing different values of v but keeping 7 and 9 r constant. The fitting 
results are shown in Fig. 3(a) with solid lines. It is surprising to see that the 
theoretical predictions and the experimental data are in good agreement even 
in this dynamic case where the shape of the pile is no longer independent 
of time. Moreover, in light of Fig. 0(b), the value of v chosen in Fig. 0( a ) 
is found to be an increasing function of H, which is qualitatively consistent 
with our experimental observation. In fact, the corresponding value of kinetic 
energy v 2 /2, as expected by conservation of mechanical energy, is decreasing 
linearly with a decrease in H for the relative larger H cases and may approach 
a constant as H — > 0. This result implies that the value of v is dominated by 
the external disturbances from the falling grains through the distance H if 
the latter is larger enough; otherwise, it becomes essentially independent of 
H, and thereby further justifies our choosing different values of v for different 
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growth stages. 

In summary, we have solved the BCRE equations analytically with a 
simple transformation. The analytical solution has been used to describe 
the surface shape of the granular pile successfully in both the static and the 
dynamic cases. This expression for the surface shape not only predicts a 
logarithmic tail as in the early work pQ, but also reveals a clear physical 
picture for its formation. 
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Granular 
materials 


Shape 


Angle of 
repose 6 C 


size 
a (mm) 


Existence of 
curved tail 


Condition of 
pile surface 


Glass 


Spherical 


30° 


= 0.3 


No 


Smooth 


Glass 


Non-Spherical 


39° 


1.25 


Yes 


Rough 


PMMA 


Spherical 


37° 


0.45 ~ 0.5 


Yes 


Smooth 


Sugar 


Cubic 


45° 


= 0.8 


Yes 


Rough 


Sand 


Irregular 


39° 


= 0.3 


Yes 


Smooth 



Table 1: Summary of the experimental results relative to the static granular 
piles. Note that the typical size of the sugar grains is obtained by measuring 
the volume of the cubic grains and calculating the typical size as the cubic 
root averaged over 10 different grains. The values of the angle of repose 
are obtained by pouring the grains in the cell mentioned earlier and then 
measuring the resulting angles of piles. The cell gap is smaller (0.25cm) than 
that used in Ref. jSj so that the values of the angle of repose here are greater 
than those in Ref. [3] due to the extra friction induced from the vertical walls. 
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Figure captions 

Fig 1 Pictures of static granular piles, (a) Spherical glass bead pile; (b) 
Irregular sand grain pile; (c) Cubic sugar grain pile. 

Fig 2 The surface shape given by Eq. (|14j) is in good agreement with 
the experimental result using sugar grains with 6 C = 45°, v = 11.5cm/s and 
7 = 23/s. 

Fig 3 (a) The surface profiles of the sandpile at different growth stages. 
The scatters are the experimental data and the solid lines represent the 
theoretical results given by Eq. (|T4*Jl with 6 C = 39°, 7 = 20/s and different 
values of v. Note that t — represents the start point of pouring the sand 
into the cell, (b) The value of v 2 /2 with v being the velocity of the gains 
within the rolling layer used to fit the pile profiles at different growth stages is 
decreasing with a decrease in H, the corresponding distance from the outlet 
of the funnel to the top of the sandpile. 
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Figure 1: 
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